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PROCESS FOR RECONSTRUCTION OF A TOMOGRAPHIC IMAGE BY AN 
ANALYTICAL METHOD COMPRISING IMPROVED MODELLING OF THE 

MOVEMENT OF THE OBJECT 

DESCRIPTION 

The object of this invention is a process for 
reconstruction of a tomographic image, of the 
analytical type, and in which improved modelling of the 
movement of the object is resorted to so as to reduce 
the artefacts of the image. 

The reconstruction of images in tomography 
implies the utilisation of radiation passing through 
the object. In numerous cases, the radiation is 
irradiation radiation partially attenuated by the 
points of the object it passes through; it also happens 
that the radiation lines are artificial and correspond 
simply to collimation lines of the detectors, which 
register an emission of particles in the object. The 
first of these processes are radiography processes, and 
the others are emissive processes, reserved more for 
living beings and in which the radiation is produced by 
an emissive body previously absorbed. Although the two 
families of processes are completely different, the 
majority of the processes for reconstruction of images 
apply to both, and this is also the case with the 
invention. 

Divergent radiation is often preferred for more 
easily enclosing the object studied.. In the case of 
irradiation radiation attenuated by the object, this is 
achieved by having a source emissive at a punctiform 
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focal point near the object and a network of detectors 
on the opposite side of the object, all of which are 
collimated towards the source; in the case of radiation 
emitted by the object, a similar network of detectors 

5 is collimated to a un focal point which corresponds 
geometrically to the preceding punctiform source, but 
to any material object. In all cases, the focal point 
and the network of detectors are moved about the object 
by taking successive views thereof, and for each of the 

10 views, the detectors measure totals of the attenuation 
or emission property along collimation lines, known as 
projections of the image of the object. When a 
sufficient number of views has been taken, there is a 
large number of projections criss-crossing through the 

15 object. The mathematical problem known as system 
inversion provides the attenuation or emission property 
which serves to form the image at each of the points of 
the object from sums of properties on the projections. 
It can be represented by a linear system of equations 

20 where the values known at the outset are the 
measurements taken by the detectors at different 
locations and the unknown ones are the values of the 
property at different points of the object. Certain 
methods of reconstruction known as algebraic 

25 effectively resort to inversion of this system; yet 
there are other methods known as analytical, where the 
value of the property at each point is calculated 
directly from . a mathematical combination of the 
projections. The present invention is one such. 

30 An important domain of tomography is the study 

of living beings for medicine. As seveiral successive 
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views must be taken, the object in question is capable 
of moving between them according to a general law of 
evolution. The different views thus represent different 
states of the object, and artefacts of reconstruction 
5 can appear only on the reconstructed image . 

Processes for avoiding this permanent problem 
consist of taking all the views at once with as many- 
sources and networks of detectors, or taking views of 
the object only in states identical to the latter, 
10 which is possible if its movement is periodical as are 
certain physiological activities such as heartbeat or 
breathing; the first of these processes is however 
expensive and the second is delicate to implement 
conveniently . 

15 This is also why numerical models of the 

movement of the object reflecting the general law of 
evolution of the object were incorporated in the 
process of reconstruction. The model generally goes 
back to a collection of shift fields of the points of 

20 the object at respective times where the views are 
taken. This can be established by other measures or by 
certain hypotheses. Several genres of models have 
already been put forward, but they are generally 
insufficient. US patent A 5 287 276, applied to the 

25 study of the thoracic cage, considered translations and 
dilations of its content, but this model is 
insufficient for other movements, such as that of the 
heart, which comprises torsion. In addition, the sets 
of projections utilised to apply the inversion formula 

30 are not very well known, which can allow new artefacts. 
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The process of the invention offers the 
advantage of using a more elaborate but also very- 
simple model of the movement of the object, comprising 
new classes of movements, and returning to a new and 
5 perfectly exact inversion formula . 

In addition, it is easy to obtain the sets of 
projections necessary and only necessary for 
reconstruction of the points of the image. 

Generalisations to more complex situations are 
10 also proposed. 

In its general form the invention concerns a 
process for reconstruction of a tomographic image of an 
especially mobile and deformable object, the image 
being a set of values of a property taken by points of 
15 the object, comprising the use of: divergent radiation 
from a focal point and passing through the object, the 
focal point being mobile about the object; an 
analytical model of mobility and deformation of the 
object defined for each position of the focal point; 
20 and an analytical calculation process for obtaining 
said values from totals of the values of the property 
along projection lines leading to the focal point and 
passing respectively by the points; characterised in 
that the model is improved, and is a variable 
25 combination being acquired, this combination comprising 
translations, rotations and homotheties of the object 
from an origin, and in that the process of analytical 
calculation comprises the following stages: 

- weighting of the measurements, this weighting being 
30 dependent on the analytical model of mobility and 

deformation of the object; 
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- derivation of the measurements weighted following 
the trajectory of the focal point considering a 
direction adapted to the model, this direction being 
kept constant, and obtaining modified measurements; 
5 - retroproj ection of the modified measurements. 

Prior to weighting, filtering of the 
measurements acquired by a Hilbert filter is often 
applied, since it is usual in other processes. The 
weighting and derivation then apply to the filtered 
10 measurements. 

The following figures are introduced to explain 
the invent ion : 

Figure 1 illustrates taking the measurements of the 
object in question, 
15 Figure 2 illustrates an equivalent diagram, the object 
being maintained in a state of reference. 
Figures 3a, 3b and 3c illustrate a complex case of 

deformation of the object. 
Figure 4 is similar to Figure 2 in the complex case, 
20 Figure 5 illustrates taking measurements with conical 
radiation, and 
Figures 6 , 7 and 8 are organigrams of three 
embodiments of the invention. 
Reference is now made to the figures. The first 
25 of them illustrates the classic situation in tomography 
of a source S of radiation and a tuned mobile network 
of detectors D on circular trajectories concentric to 
opposite positions around an object E to be studied. 
Radii R join the source S to the respective detectors 
30 of the network D. Under consideration here is a problem 
of geometric plan, which corresponds to the so-called 
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fan- shaped parallel measurement conditions, where study 
of the object E is made by superposed slices s. The 
processes implying conical radiation are also common 
and will be examined hereinbelow. The invention was 
5 designed to divergent radiation in general, which 
include the latter also. 

The trajectory T of the source S and opening of 
the radiation bundle are generally selected such that 
the latter includes the whole section of the object. 

10 We should consider a point P of the object E. A 

sole radius originating from the source S passes 
through it when a view is taken, and other radii R2, 
R3, etc. likewise pass through it for other positions 
S2 , S3 , etc . from the source S when other views are 

15 taken. These beams passing through the point P are 
utilised to determine the image of the point P in the 
inversion process. It is useful that the radii being 
considered are oriented in directions as varied as 
possible, so as to describe all the straight directions 

20 possible for each point of the object E. 

It should now be considered that the object E 
is moved and deformed. The point P moves inside the 
trajectory T and the radii R2 , R3 , etc. to be employed 
pass through respective positions P2 , P3 , etc. distinct 

25 from the point origin P. The problem of reconstruction 
to be resolved is identical to the artificial 
configuration of Figure 2, or the point P was taken as 
immobile reference and where everything happens as if 
the source S followed a trajectory T' having an 

30 irregular form and to the numerical definition noted by 

p5 = (and a(A)±n the real configuration of Figure 
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1) . It is obvious that this artificial diagram of the 
problem aids in its comprehension. According to the 
invention, the model of displacement and deformation 
selected for the object E is improved, composed for 
example from a variable combination over the course of 
acquisition of translations, rotations and homotheties, 
according to the formula (1) 



a2i(A)a22(>l)_ 



(1) 



where is the vectorial position of the point P 
10 relative to a reference such as the point O at the 

reference instant selected to carry out the 
calculations and reconstruction of the image, x is the 
position of the 

point 2 at another instant and especially a view-taking 
15 instant, and the coefficients a and b dependent on the 
time. A is a general parameter of the process, which 
varies between a minimal value X min and maximal value 
A max; the time is fixed by a monotone function de A, 
(optionally by A itself) , The reference instant where 
20 the image is reconstructed corresponds to A = 0. The 
attenuation or emission property to be calculated is 
noted f, and for the point P fixed by the vector x it 

is noted fM^) any A instant and /o(^o)^^ 

reference instant, where x is noted as Xq , 
25 The value of projection of the property measured 

over the entire radius R going from the source to a 
detector of the network D via the object E is given by 

the formula (2), where a ±s a unitary vector, and this 
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projection is noted as a function of this unitary- 
vector a and of the parameter A: 

Sxi^^oc)^\RdtfMi^)^ta) (2) 

5 where t is a parameter expressing the 

progression of the projection on the radius and taking 
the values included in the set of positive real numbers 
R. 

It is usual to apply a Hilbert filter to the 
10 projections in analytical methods, which is done here 
also, and the filtered projections are noted gHX 
according to the formula (3) 

Where ft is any vector, S"*" is the sphere unit in 
15 bidimensional, hn is indicated by the formula (4) 

hn is) = - \daisign{(j)e^''''^ . (4) 

To create the formulas the specialist could use 
either a classic operation of discrete convolation 
which will introduce a discrete version of the Hilbert 
20 filter, or a multiplication operation in the domain of 
Fourier, which will introduce an apodised version of 
the Fourier transform of the Hilbert filter. 

The inventors have also established the. 
following equation (5) 

25 gHM^Aln)^,^^^ (5) 

where the matrix- vector product A^^n is the 
direction orthogonal to the projection straight line at 
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the instant X passing through the mobile point 

corresponding to the point P, n ±s still any vector and 

or PoH designates parallel projections Pqh («/ s) , 
obtained by rearranging the projections effectively 
5 measured on the object in the reference state and 
filtered by the Hilbert filter according to the 
following formula (6) 

s being a parameter giving the distance at the 

10 origin. 

The inversion formula carried out is thus the following 
formula (7) 



Joy oJ Ja-v oJ ||jCo-r;i(a(;L))|p^^^ ^ 



where n* designates a unitary vector orthogonal 
15 to a radius linking the point P in question to the 
source S in the illustration of Figure 2, and where gpx 
is given in the following formula (8) 



The specialist could use a method of difference 
20 finished on two or three points to create this 
derivation step. 

This derivation is performed according to the 
parameter A of the trajectory of the focal point 
(source S) . It is specifically adapted to the improved 
25 deformation and mobility model as it applies to a 

direction A^pH by considering that the direction n is 

kept constant. So it is not the direction orthogonal to 
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an acquired radius, but the direction orthogonal to the 
equivalent radius of the artificial geometry, which is 
to be kept constant . 

A^n is the direction orthogonal to the 

projection straight at the instant A passing through 
the mobile point corresponding to point P. 

This inversion formula takes deformations of the 
object E into consideration and comprises, relative to 
other formulas, set in ordinary cases, a weighting as a 
function of the deformation of the object (by the 
determinant of AA) and of the position of the 

trajectory (by the standard between Xq and r>l(a(A)) . These 

barely constrained conditions are explained by an 
improved space of transformation, the projection 
straight lines of Figure 1 remain straight lines in the 
artificial geometry of Figure 2, such that the 
numerical problem can be resolved analytically. 
In formula (7) 



^ Fo-r;i(a(A))|| 



) 



the limits of the integral A^(Xo) designate a minimal 
set of positions X of the source S, such that on the 
object E in the reference state the directions of the 

straight lines linking Xq to r/l(a(A)) cover the entire 
interval of a semi -turn of trajectory without 
redundancy. If care is taken to thus limit the 
integral, the reconstruction formula is perfect. The 
integral of formula (7) could be classically made 
discrete by the specialist by a formula of trapezes. 



B 14364.3 JCI 



11 



10 



It is shown that the improved deformations do 
not sufficiently model all the objects which have to be 
studied in practice. Figures 3a, 3b, 3c, and 4 are 
offered as a revision of the reasoning of figures 1 and 
2. The object E evolves in complex fashion between the 
positions indicated successively in states XI, X2 and 
X3 or the position of the source S was also indicated 
along with the radius Rl, R2 or R3 leading to point P. 

In returning the object E to a reference state. 
Figure 4 shows that the radii designated R*l, R'2 and 
R'3 leading to point P and corresponding to radii Rl, 
R2 and R3 (passing through the same points of the 
object E) are no longer rectilinear. The formula (7) 



15 



20 



is thus no longer directly acceptable. However, it is 

approached by approximations, supposing improved 
deformation of the object E particular to each point P 
and valid around it. The deformation matrix 
corresponding to Ax, of the formula (1) 



a2i(>l)a22(A)_ 



L*2a). 



(1) 



25 



is thus calculated by the formula (9) 



|jr.,(Ar-('.))^r.,(^.r-fc), 
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where r^i(A,x) and F^jC^*^) compounds according to 

the main directions Xi and X2 of r(yl,3c) . Bx does not need 
to be calculated. The estimated shifts of all the 
points of the object E are thus the coefficients 
5 F(/l,x) which can concretely be introduced to shift cards 
which are read and executed at the moment of inversion. 



10 



15 



20 



25 



At each point P, it will still be necessary to 
determine the set of projections A^(jCo) which will be 

used for inversion, according to the abovementioned 
principle that an angular interval of a semi-turn must 
be covered by these projections. As the radii of the 
artificial geometry R*l, R'2 and etc, are no 

longer rectilinear, the directions of their tangent to 
the point P of intersection are considered. 

The inversion formula can thus be written 
according to the formula (10) 

' (/l,S(r-'(;o}^| 



where gpx is expressed by the formula (11) 



g,,(A,«(r-(3co),4~{ 



det^,(r-'(xo)) 



In addition, «(f ^(xo)i/l)is the direction orthogonal 
to the straight line acquired at the instant X passing 
through F"^(3cq) ; ^[fi{r~^{xQ\Xj) is a factor associated with 

deformation of the lengths on the radii, according to 
the formula (12) 



^^^^ ^"°\||^-(r-'(3c,)>i(r-(ic,),;5|J 
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Finally, L{xQ-r;f^{a(A))±s the distance from point P 
to the image of the source on the reference object, 
along the virtual radius . 

To the detriment of quality, but in the 
interests of greater simplicity, the term L(xo -F^ (a(>l)) 

can be replaced by ||^o~r;i(^('^))|| 



and the term 



by 1 



The improved model can be defined according to the 
real deformation of the object by an approximation 
10 according to an approximation criterion such as the 
criterion of lesser squares, the criterion of 
minimisation of the and standard, optionally 

complete by regularising on the gradient or the 
Laplacian. 

15 To date interest has been shown in 

reconstructions of the object E by slices, in fan- 
shaped parallel conditions. The preceding process can 
be extended with acquisitions by conical radiation 
under the conventional conditions of Figure 5, or the 

20 object E is represented in its three-dimensional 
entirety, also as for the network of detectors D, which 
is composed of a series of layers of detectors, each 
similar to that of the preceding figures; a radius R 
can be expressed by three parameters, namely the 

25 position A of the source S, an angle y which makes the 
radius R in the plane of the trajectory T relative to 
the central axis X of the bundle, and a cote q marking 
the layer of detectors in which the radius terminates. 
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The angle y replaces the direction vector a previously 
used in the formulas out of convenience . 
The following formula (13) 

5 gives the inversion which is undertaken in these 
particular geometric conditions and by making a 
Feldkamp approximation commonly used for treating the 
conical projections. The measurements are in this case 
multiplied by cos A, which is classic weighting 

10 compensating the oblique character of the acquired 
radius. This angle A is given in Figure 5 and 
illustrates the angle of the radius R with the plane of 

the trajectory T. In addition, -rA(a(/l))) is the 

distance from jc to the image of the source S on the 
15 object E in the reference state along the virtual 
radius R' , in projection on the plan of the trajectory 

T, and n{y~^{xQ\X) is the direction orthogonal to the 
straight line acquired at the instant A passing through 
r"*(x:o) on this plane. Finally, it is apparent that 

20 q{,^Xz~^{y^oi) is the side on the network of detectors D of 
the straight line, acquired at the instant A, passing 

through the mobile point corresponding to xO . 

A particular process, first described in French 
patent 01 07918 filed under number can still be applied 
25 here by way of some transformations. This process 
applies to rapid acquisitions made around the object 
by making several turns and consists of a 
reconstruction of the animated object by evaluating 
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its contents by blocks of projections taken on a 
fraction of a turn only of the source. The 
inversions of projections made on the blocks produce 
sub- images which are incorrect since they comprise 
5 only part of the measurements, but are obtained very 
rapidly and with the same rapidity create the law of 
displacement or deformation of the object by 
comparing homologous sub- images, taken for the same 
positions of spaced sources of a complete turn or a 

10 semi-turn; the sub-images are reconstructed at a 
reference instant for each of the groups of blocks 
and finally combined to give the complete image of 
the object. Here, the limits of the blocks are given 
in virtual geometry. The process is not otherwise 

15 modified. 

Finally, it is possible to improve the process 
to study an object subjected to periodical phenomena. 
Only those which are taken at phases of the phenomenon 
similar to those of the reference instant are selected 

20 as study blocks . 

Reference can be made to the organigrams of 
Figures 6, 7 and 8 to understand the invention. They 
detail three embodiments of the invention, or 
reconstruction is done respectively with compensation 

25 of movement; compensation of movement and time 
compensation; and compensation of movement, time 
compensation and consideration of periodicity of the 
evolution of the object. What has been previously 
described is applied in these processes. The voxels are 

30 evidently the points of the image of the object 
considered in the reconstructed image.. 



